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Abstract. We describe the current status of high-precision ab initio calculations of the 
spectra of molecular hydrogen ions (Hj and HD + ) and of two experiments for vibra- 
tional spectroscopy. The perspectives for a comparison between theory and experiment 
at a level of 1 ppb are considered. 



1 Introduction 

The molecular hydrogen ion (MHI) is the simplest stable molecule, containing 
just two nuclei and a single electron. Since the birth of the field of molecu- 
lar physics it has played an important role: it is on one hand an important 
benchmark system for detailed studies of energy levels [1], for collisions and 
chemical reactions between charged molecules and neutral atoms/molecules, of 
interactions with laser radiation and energetic charged particles, and for testing 
the respective theoretical descriptions. On the other hand, the MHI is also an 
astrophysically important molecule, involved in reaction chains leading to the 
production of polyatomic molecules. Over 800 publications have been written on 
this molecule in the last 35 years [2]. The large majority are theoretical studies. 

Concerning high- resolution spectroscopy of MHIs, only a limited number of 
investigations have been carried out, most of which a long time ago. Radiofre- 
quency spectroscopy of the hyperfine structure in several vibrational levels has 
been performed on trapped in a Paul trap [3] ; several low- lying fundamental 
ro-vibrational transitions of HD + have been measured using laser spectroscopy 
on an ion beam [4] , while rotational and ro-vibrational transitions of , 
and HD + close to the dissociation limit were investigated using microwave and 
laser spectroscopy, also on an ion beam [5,6]. The highest spectroscopic accura- 
cies reported so far were achieved in the experiments of Jefferts and of Wing et 
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al. [4,7], ~1-10~ 6 in relative units. Recently, also the dissociation energies have 
been obtained with accuracies ~ 1 • 10~ 6 [8]. Thus, the experimental accuracies 
have been far less than those achieved in hydrogen or helium spectroscopy. 

In the late 1990s, it was recognized that there are attractive reasons and 
many opportunities to study MHIs in novel ways and to achieve a much higher 
precision than previously possible [9]. Several techniques, not used before on 
MHIs, appeared to be applicable, including translational cooling, internal cool- 
ing, spectroscopy with reduced Doppler broadening, Dopplcr-free spectroscopy, 
high sensitivity ion detection. Novel laser systems not available at the time of 
the last precision spectroscopic studies can be used advantageously, among them 
diode lasers, quantum cascade lasers and the femtosecond frequency comb. The 
prospect of significantly improved experimental precision has also motivated us 
to develop more extended theoretical treatments of the MHI; in the course of 
these efforts, the accuracy of the energy levels has been increased by approxi- 
mately two orders of magnitude compared to previous work. 

Some of the above techniques have by now been implemented, and are re- 
ported here; the remaining appear to be feasible in the near future. These recent 
developments open up a number of novel applications of MHIs: 

(i) test advanced ab initio molecular calculations (in particular, QED contribu- 
tions) 

(ii) measure certain fundamental constants 

(iii) test concepts for the manipulation of molecules (state preparation, align- 
ment) 

(iv) sense fields (blackbody radiation [10]) 

(v) probe fundamental physics laws (e.g. Lorentz Invariance [11], time invari- 
ance of fundamental constants [12,13]) 

(vi) study electric dipole interactions between molecules [14] 

(vii) explore elastic, reactive and charge exchange collisions with neutral atoms 
and molecules at ultralow collision energies 

The successful demonstration of manipulation of MHIs at the quantum state 
level could also open up the possibility to study collisions with quantum-state 
resolution, i.e. where all parent particles are in specific quantum states. 

An attractive perspective of our work pursued under (i) is to eventually de- 
termine the ratios of electron-to-proton mass (m e /m p ), proton-to-deuteron mass 
(mp/rrid) and proton-to-triton mass (m p /mt) from a comparison between accu- 
rate experimental and theoretical energy level data. The basis for this possibility 
is the dependence of the vibrational and rotational transition frequencies on the 
fundamental constants. For fundamental vibrational and rotational transitions, 
the frequencies scale approximately as 

hv vib ~ \ZnieliiRoo , hv rot ~ {m e /n)Roo , (1) 

where fi = MiM 2 /(Mi + M 2 ) is the reduced mass of the two nuclei and 

is the Rydberg energy The precise dependencies have been computed in refs. 



Spectroscopy of the Hydrogen Molecular Ions. 



3 



[12,15,16]. The mass ratios m p /md 7 m p /m t and m e /m p are conventionally deter- 
mined by Penning ion trap mass spectrometry on single particles or by electron 
spin resonance of single hydrogen-like ions in a Penning ion trap. Relative accu- 
racies are currently 2.0 • HT 10 [17], 2 • 1CT 9 , and 4.6 • 1CT 10 [17,18], respectively. 
Note that in case of m e /m p , the determination involves the use of QED [19]. 
Clearly, the corresponding accuracies of v v n,, v rot represent the goal levels for 
our ongoing experimental and theoretical efforts on and HD+. 

Several aspects support the expectation that such accuracies can be reached 
in the near future. First, the lifetimes of vibrational levels are long, the shortest 
ones occurring for low-lying levels in HD + , ~ 10 ms. The relative linewidth 
due to spontaneous decay is thus of the order or smaller than 10 -13 . Second, 
Doppler broadening can be strongly reduced or eliminated by either cooling the 
molecular ions or by performing two-photon Doppler-free spectroscopy. Finally, 
collision broadening and time-of-flight broadening can also be minimized by both 
cooling and providing a good ultra-high vacuum environment. Systematic shifts 
due to light fields, trap electric fields, and trap and environmental magnetic 
fields will need to be considered; hereby it will be helpful that these influences 
be calculated accurately, using the relative simplicity of the MHI. The theoretical 
determination of the energy levels at the goal accuracy level will need as input 
nuclear properties such as the proton and deuteron nuclear radii, which may be 
obtained e.g. from hydrogen spectroscopy or nuclear scattering experiments. 

In this contribution we present an overview of our theoretical and experi- 
mental results achieved recently on MHIs. Chapter 2 describes the theoretical 
approaches for a precise computation of energy levels, including hyperfine and 
QED effects and the computation of one- and two-photon spectra. Chapter 3 
presents the development and results from an experiment for trapping and spec- 
troscopy of performed at the Universite d'Evry Val d'Essonne. Chapter 4 
summarizes an experiment on HD + at the University of Diisseldorf. 

2 Ab initio theory of H+ and HD+ 

The dissociation energies of 462 states in and 619 in HD+ in a wide range of v 
and L, vibrational and rotational quantum numbers, have been calculated some 
time ago by R.E. Moss [20,21] with a relative accuracy of ^5-10~ 9 (including 
the leading order relativistic and radiative corrections). Later the numerical 
precision of the nonrelativistic energies have been improved up to 10~ 15 — 10~ 24 
a.u. [12,15,22-26] by using variational methods. The ultimate accuracy of <~ 
10~ 24 a.u. has been obtained for the FPj" ground state [26]. These calculations 
demonstrate that at least the nonrelativistic ro- vibrational transition frequencies 
can be determined with an uncertainty well below the 1 kHz level. In this section 
we describe the calculation of QED corrections as an expansion in terms of a, the 
fine structure constant. The numerical method exploits a variational approach 
based on the Slater-type exponents as basis functions. We demonstrate that 
the frequencies of ro-vibrational transitions can be obtained in this way with a 
precision better than 1 part per billion (ppb). 
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2.1 Variational expansion 

The variational bound state wave functions are calculated by solving the three- 
body Schrodingcr equation with Coulomb interaction using the variational ap- 
proach based on the exponential expansion with randomly chosen exponents. 
Details and the particular strategy of choice of the variational nonlinear param- 
eters and basis structure that have been adopted in the present work can be 
found in [23]. 

Briefly, the wave function for a state with a total orbital angular momentum 
L and of a total spatial parity n = (— 1) L is expanded as follows: 

^ M (R,n)= Y, y l L ll A}(Ah)Gfj 2 (R,r u r 2 ), 

h+h=L 
N 

Gti 2 (R,n,r 2 ) = ]T [c n Rc[e- a " R -^-^] ( 2 ) 

71=1 

+D n l m [ e - tt - R -^ ri -^ r2 ] }. 

Here J^aKR, = R ll r l i{Yi 1 ® Y( 2 }lm are the solid bipolar harmonics, R is 
the position vector of nucleus 2 relative to nucleus 1, and ri, r 2 are positions of 
an electron relative to nuclei 1 and 2, respectively. The complex exponents, a, 
/3, 7, are generated in a pseudorandom way. 

When the exponents a n , /3„, and 7„ are real, the method reveals slow con- 
vergence for molecular type Coulomb systems. The use of complex exponents 
allows to reproduce the oscillatory behaviour of the vibrational part of the wave 
function and to improve convergence [27,23]. 

The advantage of choice (2) is the simplicity of the basis functions. It al- 
lows evaluating analytically matrix elements of the Breit-Pauli Hamiltonian and 
the leading-order radiative corrections and, more importantly, to treat in a sys- 
tematic way the singular integrations encountered in higher-order contributions 
[28]. 

2.2 Leading-order relativistic and radiative corrections 

Relativistic corrections of the leading Rooa 2 order, the Breit-Pauli Hamiltonian, 
are well known and can be found in many textbooks [29,30]. The nuclear fi- 
nite size effects are considered as contributions to this order. Details, relevant 
particularly to the case of the MHIs, can be found in [31]. 

In what follows we assume that the nuclear charges are Z\ = Z 2 = Z=\ and 
nuclear masses are denoted by capital M. The units adopted arc (h = e = m e = 1 ) . 

The radiative corrections of an order R^a 3 for a one electron molecular 
system can be expressed by the following set of equations (see Refs. [32-34]). 
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<[J,[#o, J]]/2) 

is the Bethe logarithm. The latter quantity presents the most difficult numerical 
problem in computation of QED corrections for the three-body bound states. In 
[34,35] the calculations for a wide range of ro- vibrational states in and HD + 
have been performed to an accuracy of about 7 significant digits. The operator 
J in (4) is the electric current density operator of the system 1 . The last term, 
Q(r), in Eq. (3) is the mean value of a regularized operator introduced by Araki 
and Sucher [36] for the l/(47rr 3 ) potential: 



0(r)= i im o(%^ + ( 1 ^ + 7^(r) 



(5) 



The values of this matrix element for ro- vibrational states are calculated in [31]. 

The remaining contributions in this order can be obtained from the Pauli 
form factor of an electron (anomalous magnetic moment): 



£ (3) = na 2 Z 



(?) 



1 /a 
2 



<5(r 1 )+5(r 2 )}. 



and from the one-loop vacuum polarization: 

4a 3 Z 



'-'vp 



(<J(n)+<J(r 2 )) 



(6) 



(7) 



2.3 RacCx 4 order corrections in the nonrecoil limit 

The contribution of recoil corrections, proportional to (m/M), in the Rocoft 
order are too small for our present consideration and may be neglected. Radia- 
tive corrections for a bound electron in an external field are known analytically 

1 J~^2 a z a p a /m a , where z a , p a , m a are the charge, impulse, and mass of a particle 
a. The sum is performed over all particles of the system. 
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Fig. 1. Adiabatic potential of the ma 6 order contribution to the Dirac energy of the 
two-center problem (Zi = Z 2 = 1). Dashed curves are the squared densities of the 
ground and first vibrational state wave functions of H J ion. 
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The last equation includes both the Dirac form factor and polarization operator 
contributions. 

The i?ooa 4 relativistic correction is obtained using the adiabatic "effective" 
potential for an ma e term in the a expansion of the two-center Dirac energy 
(see Fig. 1). Averaging over the squared wave function density of a state one 
gets E$. The adiabatic potentials have been obtained recently with about 5 
significant digits [37], and the Born-Oppcnhcimcr solution yields approximate 
wave functions at the (m/M) w 10~ 4 level, which allows to claim that E$ is 
now known to 4 digits. 

Some higher order radiative corrections for a bound electron in an exter- 
nal field are also known in an analytic form [38,39] and can be included into 
consideration: 



Eg =a"Z 3 \n\Za)- 2 [-l] {5^ + 5^)). 



(9) 



The electron ground state wave function may be approximated by tp e (r e ) = 
C[ipi s (ri) + ^is( r 2)], where ipi s is the hydrogen ground state wave function 
and C is a normalization coefficient. Thus, one may use this approximation to 
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evaluate other contributions in the R^a 5 order: 



A 61 \n(Za)- 2 + A 60 (5(n)+5(r 2 )> 



E. 



.(B) 
2loop 



a" 

7T 



Z 2 [S 50 ] (J(n)+<f(r 2 )) 



(10) 



where the constants A6i, Ago, and £50 are taken equal to the constants of the 
Is state of the hydrogen atom A 61 = 5.419... [40], A 60 = -30.924... [41], 
and -B50 = —21.556 . . . [42] (see also Ref. [39] and references therein). The final 
theoretical uncertainty in the transition frequency (see Table 1) is determined 
by the total contribution of the last two equations. 



2.4 Hyperfine structure of states 

The leading order contribution to the hyperfine splitting of the ro-vibrational 
states is calculated using the spin-dependent part of the Breit-Pauli interaction 
Hamiltonian, with phenomenological values for the nuclear magnetic moments 
and the electron anomalous magnetic moment. The hyperfine levels of HD + , 
E-vlfsj, are labelled with the quantum numbers F, S and J of the intermediate 
angular momenta F = I p +s e , S = F+I^ and of the total angular momentum J = 
L + S [43] . In case of due to Pauli exclusion principle the total nuclear spin I 
is uniquely defined by L and parity of the electronic state. The following coupling 
scheme is adopted: F = I + s e and J = L + F [44]. The hyperfine structure (HFS) 
of the ro-vibrational states of HD + consists of 4, 10 or 12 hyperfine sub-levels 
for L = 0, L = l and L>2, respectively (see Fig. 2). The multiplicity of the HFS 
of is reduced to 1 for L = 0, 5 for L = l, 2 for even, and 6 for odd L states. 
Typically, the hyperfine splitting of the lower ro-vibrational states of HD+ and 

is about 1 GHz. The uncertainty in the hyperfine spectrum is related to the 
unknown contribution of the spin interaction terms of orders 0(R OD a 4 (m/M)) 
and higher, which have not yet been taken into consideration, and is estimated 
not to exceed 100 kHz. 

Each transition line between ro-vibrational states is split into a multiplet of 
hyperfine components, corresponding to the allowed transitions i — ► / between 



Table 1. Summary of contributions to the (v = 0, L = 0) 
frequency (in MHz). 



■ (v' = 1, L' = 0) transition 



AEnr 

AE a2 
AE a3 

AE>, 



HD+ 



65 687 511.0686 
1091.041(03) 
-276.544(02) 
-1.997 
0.120(23) 



57 349 439.9717 

958.152(03) 
-242.118(02) 
-1.748 
0.106(19) 



AEtot 65 688 323.688(25) 57 350154.368(21) 
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Fig. 2. Hyperfine structure of a ro-vibrational state of HD + with L > 2. 



the states of the hyperfine structure of the initial and final states. Whether these 
hyperfine components will be resolved or the spectrum profile depends on the 
initial and final state lifetime and on the experimental conditions (transition 
linewidth ry, laser intensity /, temperature, interaction time etc.) Examples of 
spectral data are presented in experimental sections of our review (see Fig. 13). 
The shape of the profile also depends on the population of the initial hyperfine 
states. The general expressions for the probability per unit time for one- and 
two-photon transitions between ro-vibrational states of the MHIs with account 
of the hyperfine structure are given in Ref. [45]. 

The probability per unit time for the hyperfine transition i — > / at resonance 
(averaged over the magnetic numbers of initial and final states), 1/^, may be 
represented in the form: 



r r , r r 1 (v'L'\\d\\vL)* | f 



Here T v iu,vL is the probability per unit time of laser-stimulated dipole transi- 
tions between ro-vibrational states, d \ \vL) is the reduced matrix element 
of the electric dipole moment of the HD + ion d = ^ a z a r a , and 



!>,< = V(2J' + 1)(2L + 1) £ (-lf'+ J + L ' { L I, L j } (3l„ s „ 

F"S" ^ J 



&»s», (12) 



where p^ffi are constant amplitudes of the state vectors of the hyperfine 
states: 

\vLFSJ, J z ) = P^ F S" J E cJ lm,S»<*vLm{R, r!) X (F"5", C), (13) 

F"S" M( 

determined from the effective Hamiltonian of spin interaction. Here £) 
are basis spinors of definite values of F, S and S z in the space of the spin 
variables. The relative intensity of the hyperfine components of a transition line 
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Fig. 3. Hyperfine splitting and intensities of the two-photon ro-vibrational transition 
line for the H+ ion: (a) (v = 0, L= 1) -> (v'= 1, L' = 1), (b) (v = 0, L = 2) -> (v' = 1, L' = 2). 

between ro-vibrational states is thus determined by the amplitudes Tfj. In case 
the individual hyperfine components cannot be resolved, the observable intensity 
is reduced to the intensity of the dipole ro-vibrational transition A'L'.tjL, in 
agreement with the identity J2f ^/,» = 1- 

The hyperfine structure of the one- and two-photon transition lines includes 
a large number of components, most of which, however, are suppressed. There 
are as well dominant (or "favoured') transitions between states with similar 
spin structure, such as (vLFJ) — ► (v'L'FJ') with A J = AL (for H^~). In such 
pairs of homologous hyperfine states the spin-dependent corrections to the ro- 
vibrational energies E v > l> fji and E v lfj have close values, which partially cancel 
each other when evaluating the spin correction to the resonance transition fre- 
quency (Ev'l'fji —E V LFj)/h. Indeed, the favoured hyperfine transitions span 
over a frequency interval less than 25 MHz (see Fig 3). It is natural to expect 
that the unknown contributions to the frequency of the favoured transitions from 
the spin interactions of order i?ooa 4 (m/M) and higher also tend to cancel each 
other; therefore, the theoretical uncertainty of the resonance frequency of the 
favoured hyperfine sublines will be less than ~ 5 kHz. 

2.5 Two-photon transition probabilities 

In order to assess the feasibility of Doppler-free two-photon spectroscopy in H^" 
or HD + , it is essential to evaluate transition probabilities between ro-vibrational 
states. This can be done using the formula from second-order perturbation the- 
ory [46], and the accurate wavefunctions provided by variational calculations. 
Assuming that magnetic sublevels are equally populated, the two-photon tran- 
sition probability at resonance between states v,L and v',L' is 




(14) 



where I is the excitation intensity, // the transition linewidth, and 




(15) 
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Q(°\ Q( 2 ) are respectively the scalar and tensor parts of the two-photon transi- 
tion operator 

Q = de ir-u de ( 16 ) 

47re ag E - H 

Here, E = [E(v, L) + E(v' ', L')}/2 is the one-photon resonance energy and e the 
exciting field polarization. The two-photon transition probabilities were calcu- 
lated in [22,47]. 

For the case [47] there exists a quasi-selection rule Av = ±1, and the 
dimcnsionlcss transition probabilities Q v L v , L , are rather small, of the order 
of 1, which is due to the level structure of H^~. If we consider the example 
of transitions between L = states, these states are of 1 S e symmetry, and 
there is no resonant intermediate level of 1 P° symmetry that could enhance the 
transition probability. 

The situation is different in the HD + case [22], since there is no splitting 
between singlet and triplet state due to the loss of exchange symmetry between 
the nuclei. As a result, for a transition between L = (S e ) states, there exist 
intermediate bound P° levels which can be very close in energy and efficiently 
enhance the transition probability. This is most likely to happen if the difference 
between v and v' (the initial and final vibrational quantum numbers) is an even 
number. In this case, the state (v" = (v + v')/2,L" = L±l) is often close to 
the middle energy E = [E(v,L) + E(v' , L')]/2, because of the quasi-harmonic 
structure of vibrational levels. As a result, some of the most intense two-photon 
lines arc Av = 2 transitions in the 5-6 /xm range or Av = 4 transitions in the 2.5- 
3 (im range, accessible e.g. with continuous-wave optical parametric oscillator 
or quantum cascade lasers. The dimensionless transition probabilities Q v L v i L i 
can reach values as high as 300 for the (v = 0, L = 1) — > (v = 2,L = 1) transition 
at 5.366 ^m. Thus, the HD + molecular ion is a promising candidate for precise 
two-photon spectroscopy. 

The next step is to consider the hyperfine structure of two-photon transition 
lines. The representation of the hyperfine state vectors in the case is 

\vLFJJ z ) = £ f3^ FJ £ C J L J ^ F , C Vvlm(R, ri) C). (17) 

F> MQ 

The definitions for x(i 7 ",C) an d (3jfi FJ are similar to those in (13). Assuming 
equal populations for hyperfine magnetic sublevels, the dimensionless transition 
probability between levels \i) = \vLFJ) and [/) = \v'L'F'J') is: 



Q iJ = (2J'+l)J2 



(^||Q( fe )||^E F "(-l) J ' +IH - F ''{j,^j}/3^^" 



2k + 1 

fe=0,2 



(18) 



The hyperfine structure of two different transitions for linear polarization is 
shown in Fig. 3. Transitions between odd L states (Fig. 3(a)) comprise between 
25 to 34 components, 5 or 6 of which are favoured. Transitions between even 
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Fig. 4. Born-Oppenheimer electronic energies (in a.u.) of the adiabatic potential terms. 

L states (Fig. 3(b)) have a much simpler structure, the total nuclear spin be- 
ing zero. The non zero even L spectrum comprises only two main components 
(verifying A3 — 0) together with two weak satellites. Transitions between L = 
states are structureless, which makes them especially attractive from a metro- 
logical point of view. 



3 Two-photon spectroscopy of H+ 

A two-photon vibrational spectroscopy experiment aimed at the determination 
of the electron to proton mass ratio is being setup at the Kastler Brossel Labo- 
ratory. We begin by recalling the basic spectral features of the MHI, and discuss 
the planned experimental sequence. In the second part, we report on the present 
status of the experimental setup. It is composed of a hyperbolic Paul trap in 
which a few thousand ions can be confined, a UV laser for ion preparation 
and detection by state-selective photodissociation, and a narrow-line, tunable 
laser system that will excite the two-photon transition. 

3.1 H2 level structure 

Although the Born Oppcnhcimcr (BO) approximation is not relevant for highly 
accurate calculations, it remains a very convenient tool to get a useful insight 
into the level structure. In order to understand the processes discussed here, 
it is enough to consider the first two BO electronic curves: the ground state lsa g 
and first excited state 2pa Ul which are depicted in Fig. 4. 

The exact symmetries of the system are the total spatial parity ir and the 
exchange of nuclei P12; the g/u electronic parity 7r e used in the BO approach is re- 
lated to them by 7r e = 7rPi 2 . The bound levels of can be labelled v, ( 2I+1 L e '°) 
where v and L are the vibrational and orbital quantum numbers, I is the total 
nuclear spin quantum number and (e,o) stands for the total parity. Since the 
total spatial parity is n = (— 1) L , the lsa g curve only supports 1 S e , 3 P°, 1 D e . . . 
levels. 



12 Bernhard Roth et al. 



0,01 
0,0001 

0,01 
0,0001 

s « 0,01 

2 o.oooi 

D 

0,01 
0,0001 

50 100 150 200 250 300 

X (nm) 

Fig. 5. Photodissociation cross sections of the L = 0, v levels of H^. The dashed lines 
are the result of a Born Oppenheimer calculation [50], which takes the lsa g and 2pa u 
electronic curves into account. The nodal structure of the cross section reflects that 
of the vibrational wavefunction. The solid lines are obtained from an exact variational 
method which fully takes into account the three body dynamics [51]. The additional 
nodal structure appearing in the short wavelength domain can be interpreted as the 
photodissociation to higher excited electronic states (3da g , 2pir u , 4fa u , . . .). The dot- 
dashed vertical line corresponds to the KrF laser wavelength of 248 nm. 

The 2pa u electronic curve presents at large internuclear distances a weak 
attractive potential that supports two bound L = energy levels [20,48]. Some of 
those states have been observed by microwave or laser spectroscopy experiments 
[49]. At higher L the 2pa u potential supports bound or dissociative 1 P°, 3 D e 
1 F° . . . states that can be calculated numerically using either the variational, or 
the complex coordinate rotation method. 

3.2 One-photon transitions: photodissociation 

The selection rules for one-photon dipole transitions are AL = ±1 and AI = 0. 
As a consequence, transitions between bound ro- vibrational states of are for- 
bidden (in contrast with the HD + case), resulting in very long-lived states. On 
the other hand, one-photon photodissociation transitions from lsa g to 2pa u elec- 
tronic states are allowed. The photodissociation cross sections a v of the ( 1 S e ,v) 
states have been first computed by Dunn [50] in the Born-Oppenhcimcr approx- 
imation, and then using the perimetric coordinate variational method in [51]. 
The results are given in Fig. 5. They show that a laser source in the 250 nm 
range can be used to selectively photodissociate the v — 1,2,3,... vibrational 
states while keeping the ions in the v = level since the successive cross section 
ratios a v+ i/a v are 214, 40, 10, for v = 0, 1, 2, respectively. 

3.3 Two-photon transitions: choice of the transition 

One-photon transitions between bound states being forbidden, a high resolution 
study of the vibrational structure of is only possible using Doppler-free two- 
photon spectroscopy. Two-photon transitions obey the selection rule AL = 0, ±2 
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as well as the quasi selection rule Av = ±1, as discussed in section 2.5. The 
corresponding (v,L) — > (v' = v + l, L') transition frequencies lie in the 8-12 /j,m 
range. Among them, we have chosen to probe the v = — > v = 1 transitions, for 
L and L' equal to or 2 and eventually 1 or 3. We now give the arguments that 
explain this choice. 

The first condition to fulfill is that it should be possible to prepare a large 
enough number of in the initial state of the transition. The ro-vibrational 
populations of ITj~ ions, after creation by electron impact ionization of a low pres- 
sure H 2 gas, have been studied both theoretically and experimentally [52]. The 
vibrational populations are linked to the overlap of the H 2 and vibrational 
wave functions (Franck-Condon principle); they are found to be of the order of 
12, 18, 18, 15, 11, 8, 5, 4% for the first few levels. The rotational populations of 
the ions are those of the H 2 mother molecules, i.e. 12, 28, 28, 18, 8% at 300 
K. Moreover, we have shown in the previous paragraph that UV photodissocia- 
tion provides a convenient way to prepare ions in the ground vibrational state; 
it is then desirable to choose a v = state as initial state of the transition, with 
L between and 3, L= 1 or 2 being the best choices with respect to the number 
of ions. The same photodissociation process can be used to detect the ions in 
the excited v = l state. 

The hyperfine structure of two-photon lines should also be considered; it is 
apparent from Fig. 3 that it is simpler for transitions between even L states. 
Interpretation of experimental data is likely to be easier in this case. 

The intensities of the various (v = 0,L) — ► (V = 1,L') (with low L) two- 
photon lines are of the same magnitude; the choice of a particular transition 
depends mostly on the availability and characteristics of laser sources at the 
required wavelength. The whole mid-infrared range is accessible by the recently 
developed quantum cascade lasers (QCL); L' — L transitions are especially at- 
tractive, because they lie within the spectrum of C0 2 lasers (A ~ 9-10 /im). 
Also, a number of frequency reference molecular absorption lines are known in 
this range [53] . The first transition that is going to be probed in our experiment 
is the (v = 0, L = 2) — -> (v 1 = 1, V — 2) line at 9.166 nm. The details of coinci- 
dences with C0 2 lines and molecular reference lines, which make this transition 
favorable, are explained below. 

3.4 Experimental sequence 

The two-photon transition matrix elements \Q v ,l,v'.l' 2 of the "favoured" hyper- 
fine components of two-photon transitions are of the order of 0.2 (see Fig. 3). A 
typical QCL can deliver about 50 mW of single-mode optical power. Assuming a 
perfect coupling to a build-up cavity of finesse 1000 with a waist of 1 mm, one ob- 
tains a laser flux of 15 W/mm 2 . Assuming an instrumental width jy wlO kHz, 
equation (14) yields transition rates of about 70/s. This order of magnitude 
shows that long interaction times are needed and that one has to work with a 
cloud of trapped ions, having a radius of the order of the beam waist. 

The considerations of the previous paragraph show that vibrational two- 
photon spectroscopy of can be performed by (2 + 1') resonance enhanced 
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Fig. 6. (a) Simulation of the experimental setup for ion creation, trapping and detec- 
tion using the SIMION7 software. MCP is a multichannel plate detector, (b) Time of 
flight spectrum showing the H + , and H3 species confined in the Paul trap. 

multiphoton dissociation (REMPD). This process is very similar to that im- 
plemented for HD + vibrational spectroscopy and described in more detail in 
Section 4. 

The experiment will be conducted in the following stages: 

• simultaneous creation, trapping and selection of (L, v = 0) H^~ ions. 

• excitation of the (L,v = 0) — > (L',v' = l) two-photon transition. 

• photodissociation of the (L,v=l) Hj ions. 

• time of flight detection of H + and H^" ions. 

3.5 Experimental setup 

The ion trap is depicted in Fig. 6(a). It is a hyperbolic Paul trap with a ring of 
inner radius ro = 8.2 mm and two end caps separated by 2z = 6 mm. Two pairs 
of holes (2 mm in diameter) arc drilled in the ring along two orthogonal axes to 
shine the ion cloud with the UV and IR light. Both end cap electrodes are AC 
grounded. A RF voltage (about 200 V peak to peak amplitude at 10.3 MHz) and 
a continuous voltage of a few Volts are applied to the ring electrode, resulting 
in trapping well depths of a few eV. 

The ions are produced by electron impact ionisation from the residual 
H2 gas. The electron gun is made of a tungsten wire and a Wehnelt cylinder; 
it is typically turned on for 100-200 ms. A 1 mm hole in one of the end cap 
electrodes allows access to the trap. 

The contents of the trap are analyzed by applying a short negative high 
voltage pulse to the second end cap, thus extracting the ions from the trap 
through a 2 mm hole. The extracted ions are accelerated and focused onto a 
multi-channel plate (MCP) detector located 7 cm away, a long enough distance 
to separate by time of flight the H + , H^ and H^ ions that are simultaneously 
produced and trapped. A typical time of flight spectrum is shown in Fig. 6(b). 
Up to a few thousand Hj ions can be stored in the trap. The ion lifetime is of 
a few seconds and is limited by the residual pressure in the vacuum chamber. 

The undesirable H + and Hjj~ ions arc eliminated using the parametric exci- 
tation of their secular motion, by applying RF voltage in the MHz range on one 
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Fig. 7. (a) Absorption spectrum of formic acid (HCOOH) [53] . The line intensities 
are given in cm -1 /(molecule cm -2 ), (b) Two-photon transition probabilities in atomic 
units. The central peak is made of two close components (see Fig. 3). The dashed 
line is the 9R(42) CO2 emission line. The detunings between the 9R(42) CO2 line and 
the HCOOH lines are indicated in MHz. The C0 2 laser is locked to the HCOOH line 
indicated by the star. 

of the end cap electrodes during the ionisation process. A KrF excimer laser at 
248 nm is used to photodissociate the v > 1 states in order to produce a (L, v = 0) 
ions cloud. The ions are shined by 20 mJ pulses during the filling of the trap. 
The characterization of ro-vibrational populations of the resulting ion cloud is 
now in progress. 

Since all the bound states of H j" are metastable, the natural widths of the two- 
photon transitions are extremely small. In Paul traps, the ion cloud temperature 
is of the order of magnitude of the potential depth expressed in K, i.e. « 10 4 K 
in our trap. Under those conditions, the two-photon linewidth Tt (appearing in 
Eq. (14)) is expected to be limited by the second-order Doppler effect, i.e. of the 
order 10 kHz. It will limit the ultimate frequency resolution of the experiment 
at the 3-10 -10 level, and the mass ratio resolution at the 6-10~ 10 level. 

Ion cooling will thus be necessary in order to reach the metrological objective 
of the experiment at the 10 -10 level. Nevertheless, the first step of the experiment 
is the observation of a two-photon transition, which is feasible with hot ions using 
a kHz linewidth laser source. 

3.6 Two-photon excitation laser source 

The laser system we have built is aimed at exciting the (L = 2, v = 0) — > 
(L' = 2, v' = 1) two-photon transition at 9.166 /im. In this range, two kinds 
of laser sources are available. Single-mode CO2 lasers have high output power 
and sub-kHz linewidths, but are hardly tunable on ranges exceeding 100 MHz, 
i.e. much smaller than the 1.65 GHz gap between the closest CO2 emission line 
(9R(42)) and the H^ line (see Fig. 7(b)). Recently, single mode quantum cas- 
cade laser (QCL) became commercially available. They can be tuned over about 
10 cm -1 (300 GHz) through their temperature and injection current, but have 
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Fig. 8. Simplified setup of the CO2/HCOOH phase-locked quantum cascade laser 
source. The QCL is mounted in a liquid nitrogen optical cryostat. The mixer is a 
room temperature HgCdZnTe detector. Solid lines are optical paths. Dashed lines are 
electrical paths and servo loops. A.O. is an acousto-optic modulator. 



large linewidths of the order of a few MHz. Several experiments have shown 
that the lincwidth can be reduced well below the kHz level by injection-current 
locking the QCL to a molecular line [54] or to a high finesse Fabry-Perot cavity 
resonance [55]. We have developed a laser source that takes advantage of both 
the narrow linewidth of the CO2 laser and the tunability of the QCL [56]. 

The setup is shown in Fig. 8. A CO2 laser oscillating on the 9R(42) line is 
frequency shifted by 128 MHz and stabilized on the intracavity saturated ab- 
sorption signal of the (21, 3, 19) — > (21, 2, 20) line of the v§ band of formic acid 
(HCOOH) (see Fig. 7(a)). The absolute frequency of that transition (32 708 263 
980.5 kHz) has recently be determined with an uncertainty of 1 kHz [57] by sum 
frequency mixing with a 30 THz wide visible femtosecond frequency comb [58]. 
The QCL is operated in a liquid nitrogen optical cryostat. The output power 
is 50 mW with a 700 mA injection current and a temperature of 80K. The 
QCL is phase-locked to the C0 2 laser with a tunable frequency offset in the 
500-2000 MHz range [56] . The analysis of the beat note spectrum under locked 
conditions (see Fig. 9) shows that we have realized a narrow-line tunable laser 
source well suited to probe the Hj two-photon lines, and also the ro-vibrationnal 
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Fig. 9. Phase-locked beat note between the QCL and the CO2/HCOOH frequency 
reference. RBW 10 kHz, VBW 1 kHz. The spectrum shows a loop bandwidth of the 
order of 6 MHz. The central peak is extremely narrow, with a -3dB width smaller than 
the 200 Hz resolution of the spectrum analyser. The inset shows the free running beat 
note with the same scale and a 500 kHz RBW. 
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spectrum of HCOOH [58] or other molecules (NH 3 , . . . ) of atmospheric or as- 
trophysical interest. 

4 Cooling and Spectroscopy of HD+ 

In experiments performed at the University of Dusseldorf, the MHIs H^~, , 
and HD + were cooled to temperatures of ~ 10 mK in a radiofrcqucncy trap, 
by sympathetic cooling with laser-cooled beryllium ions. High-resolution spec- 
troscopic studies of several rovibrational infrared transitions in HD + were per- 
formed. Hyperfme splitting of the lines was observed, and is in good agreement 
with theoretical predictions. The transitions were detected by monitoring the de- 
crease in ion number after selective photodissociation of HD + ions in the upper 
vibrational state. 

4.1 Preparation and characterization of cold MHI ensembles 

MHIs arc just a few of a multitude of species that can be cooled to mK tempera- 
tures, by sympathetic cooling [59,60] - the molecular species and a laser-coolable 
atomic species, with the same sign of charge, are simultaneously stored in a ra- 
diofrequency trap. Laser cooling the atoms then also efficiently cools the molec- 
ular ions via the long-range Coulomb interaction. Temperatures below 20 mK 
can be reliably reached. We have shown that using Be + ions as coolant permits 
to cool sympathetically ions from mass 1 to mass 200 amu [61-63]. A heavier 
atomic coolant species can be used to extend the mass range. For example, using 
138 Ba+ as coolant molecular ions up to mass 410 amu have recently been cooled 
[64]. 

We use a linear quadrupole trap to simultaneously store both Bc+ and MHIs. 
The radiofrequency trap is driven at 14 MHz, with a peak-to-peak amplitude 
of 380 V. This results in a radial Mathicu stability parameter q r ~ 0.13 for 
HD + . The trap is enclosed in a UHV chamber kept below 10~ 10 mbar. The 
chamber is equipped with a leak valve for the controlled introduction of gases. 
An all-solid-state 313 nm laser system is used for cooling Be+ [65]. 

To load Be + ions into the trap, atoms are thermally evaporated from a beryl- 
lium wire, and ionised by an electron beam. The molecular loading is achieved by 
leaking in neutral gas at a pressure of <~ (1—3) 40~ 10 mbar, ionised by an electron 
beam with an energy of 200 eV, and a current of <~ 30 fiA, for a loading time of 
2 s. This produces mixed-species crystals like those shown in Figs. 10(a,b). The 
ions with a higher charge-to-mass ratio (in this case the molecular ions) expe- 
rience a stronger trap pseudopotential, and thus form a dark (non-fluorescing) 
core to the crystal. The asymmetric distribution of species along the z-axis ob- 
served in Fig. 10(b) is caused by the light pressure of the cooling laser on the 
beryllium ions. 

The observed crystals are reproduced by molecular dynamics (MD) simula- 
tions [61,66]. Visual matching of overall structure, structural details and blur- 
rings of CCD and simulated images allows fitting the ion numbers and tem- 
peratures of the different species. The number of ions of different species given 



18 Bernhard Roth et al. 




Fig. 10. Fluorescence images of (a) a large ion crystal with a high fraction of sympa- 
thetically cooled ions (approximately 1200 light ions and 800 Be + ions) (b) a smaller 
crystal containing approx. 690 Be + ions, and 12 (exactly) HD + ions, and simulated 
images of this crystal at (c) 20 mK, (d) 12 mK, (e) 8 mK and (f) 6 mK. In the simu- 
lations, the beryllium ions are shown in red and the HD + ions in blue. Laser cooling 
beam propagation is to the right, along the z-axis [61]. 

in Fig. 10(b) were found in this way. In the simulations we assume an ideal 
linear trap, use the quasipotential approximation and model heating effects by 
stochastic forces on the ions. The obtained temperatures are thus effective secu- 
lar temperatures. Fig. 10 shows a determination of the temperature; agreement 
is found for a Be + temperature of appox. 10 mK. This sets an upper limit, as our 
experimental images are also limited in sharpness by our detection optics, CCD 
resolution, and sensor noise, which are not considered. The temperature varies 
depending on crystal size and cooling parameters, and is typically in the range 
5 mK to 15 mK, with smaller crystals generally colder. These temperatures are 
consistent with measurements of the fluorescence lineshape of the Be + ions. 

For all species of molecular ions studied here, our MD simulations show that 
the sympathetically cooled molecular ion ensemble is also crystalline, i.e. its 
time-averaged ion distribution is strongly inhomogeneous, and that it is strongly 
thermally coupled to the Be + ions. Assuming similar heating effects for the 
molecular ions and the Be + ions, the simulations show that the moleculars ions 
have a temperature similar to that of Be + , due to the strong Coulomb coupling. 

The trapped species are identified and the time evolution of their numbers is 
monitored by excitation of their mass-dependent radial (secular) modes, using 
a spatially homogenous and temporally oscillating electric field. For HD + ions 
the measured secular frequency was «770kHz, significantly shifted from the 
calculated single-particle frequency, due to Coulomb coupling between different 
species in the trap [67]. Excitation amplitude, sweep rate, and covered frequency 
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Fig. 11. Principle of (l + l') REMPD spectroscopy of HD + ions. A tunable IR diode 
laser excites a rovibrational overtone transition (v = Q,L)— >(w' = 4, L'). The HD + ions 
excited to the v'=4 vibrational level are dissociated using cw 266 nm laser radiation: 
HD + (V = 4) + hv — > H + D + , or H + + D. Due to different Franck-Condon wavefunction 
overlap, the calculated UV absorption cross section from the «' = 4 level (~ 2.4x10" 17 
cm 2 ) is about 7 orders of magnitude larger than from v = [68]. Energy values represent 
total binding energies of the molecule [69]. 

range were chosen so that the ion crystal had sufficient time to cool back to its 
initial temperature between individual excitation cycles. The excitation heats 
both the molecular ions and the atomic coolants, which changes the scattering 
rate of 313 nm cooling light by the Be + ions. The HD + secular resonance becomes 
visible in the Be + fluorescence, and its strength is proportional to the amount 
of HD + ions in the ion crystal. 

4.2 Spectroscopy of HD+ 

The choice of HD + for spectroscopic studies was made because of the availabil- 
ity of dipole-allowed ro-vibrational transitions which simplify the spectroscopic 
techniques. Nevertheless, vibrational spectroscopy in the electronic ground state 
in near-absence of collisions, as is the case for the present molecular ions en- 
sembles, is faced with the difficulty that molecules excited to a vibrational level 
decay only slowly, implying very low fluorescence rates. As the fluorescence wave- 
lengths are in the mid to far infrared, photon counting would require a sophisti- 
cated detection system. We circumvent this difficulty by applying the technique 
of (1+1') resonance enhanced multiphoton dissociation (REMPD): the molecules 
are excited by an infrared (IR) laser and then selectively photodissociated from 
the upper vibrational state by a second, fixed-wavelength ultraviolet (UV) laser 
(Fig. 11). The remaining number of molecular ions is the quantity measured as a 
function of the frequency of the IR laser. As the molecular sample is small (typ. 
40-100 ions) the spectroscopy requires the spectra to be obtained by repeated 
molecular ion production and interrogation cycles. The lasers employed are a 
single- frequency, widely tunable diode laser at 1.4 /im (Agilent 8 1480 A), and a 
resonantly frequency-doubled Yb:YAG laser at 266 nm. The IR laser linewidth 
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Fig. 12. (a) Initial ion crystal: wllOO Be+, ^100 HD+, and ^20 D + ions at ^20 mK 
(the presence of cold HD + ions is obvious from the dark crystal core), (b) Repeated 
secular excitation of the crystal in (a) at 3 V amplitude. The excitation frequency was 
swept between 500 kHz and 1500 kHz. The IR laser is tuned to the maximum of the 
(v = 0, L = 2) — * (v' = 4, L' = 1) line. The curve is an exponential fit with a decay 
constant of 0.04 s . (c) Ion crystal after dissociation of all HD + ions: Rill 00 Be + and 
~50 D + ions at «20 mK. (d) Measurement cycle consisting of repeated probing of the 
number of HD + ions before and after exposure to the spectroscopy lasers [69] . 

was ~ 5 MHz, and its frequency was calibrated with an accuracy of 40 MHz by 
absorption spectroscopy in a water vapor cell. 

Due to the weak coupling between external and internal (rotational) degrees 
of freedom, the internal temperature of the HD + ions is close to room tem- 
perature, in thermal equilibrium with the vacuum chamber [70,10]. There is 
significant (>5%) population for rotational levels up to L = 6. Indeed, we have 
observed 12 transitions between 1391 nm and 1471 nm, from lower rotational 
levels L — to L = 6. 

The loss of HD + ions not only depends on the REMPD process, but also on 
transitions induced by blackbody radiation (BBR). We modeled the loss of HD + 
by solving the rate equations for the populations of all (v,L) levels interacting 
with the IR and UV lasers, as well as with the BBR radiation at 300 K. The 
theoretically obtained excitation spectrum (see Fig. 13 and text below) of the 
levels probed by the IR laser is included, but for the remainder of the calculation 
hyperfine structure, due to electron, nuclear and rotational spins, is ignored. 
The rovibrational transition moments involved are taken from [71]. The rate of 
dissociation by UV light is obtained using cross sections from [68] . For typical UV 
intensities, dissociation rates of 10 2 -10 3 s _1 are found. The rate equation model 
reveals two different timescales at which the HD + number declines during a 
typical experiment. A first, fast (< 1 s) decay occurs when the IR laser selectively 
promotes HD + ions from a specific (v = 0, L) level to a rotational level in t/ = 4, 
from which they are efficiently photodissociated. This process rapidly dissociates 
those (v = 0, L) HD + ions which are in the hyperfine states probed by the IR 
laser. The remaining molecular ions (a significant fraction of the total initial 
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number) are dissociated significantly slower, essentially at the rate at which 
the hypcrfinc levels of (v = 0, L) are repopulated by BBR and spontaneous 
emission. For example, for the (y = 0, L = 2) — > (v' = 4, L' = 1) transition, 
and for typical intensities of 6 W/cm 2 for the IR and 10 W/cm 2 for the UV 
laser, the fast HD + decay takes place at a rate ~10 s _1 (which is not resolved 
experimentally), whereas the decay due to BBR-induced repopulation occurs at 
a rate of ^0.04 s _1 . The latter rate is fairly consistent with the measured decay 
depicted in Fig. 12(b), but observed decay rates depend strongly on which part 
of the hyperfine spectrum is interrogated. This points at a shortcoming of the 
simple rate equation model used here, and our observations can probably be 
explained precisely only by a rate equation model which takes the full hyperfine 
structure of all involved (v, L) levels into account. 

As an example, Fig. 12(b) shows the time evolution of the HD+ secular 
excitation resonance while the HD + ions are excited on the maximum of the 
rovibrational line (v = 0,L = 2) -> (v' = 4,1/ = 1) at 1430.3883 nm. The de- 
crease of the HD + resonance in the secular excitation spectrum, induced by the 
REMPD process, is accompanied by a decrease of the dark crystal core con- 
taining the MHIs. The secular excitation spectrum also shows an increase of the 
number of D + ions, which result from the dissociation of excited HD+ ions. These 
ions are sympathetically cooled and remain in the crystal core. Fig. 12(c) shows 
the mixed-species ion crystal after all HD+ was dissociated. The dark crystal 
core has shrunk significantly, and the crystal now contains «1100 Be + and «50 
D + ions. Assuming equal probability for photodissociation to D + and H + , this 
number indicates that most generated D + ions are sympathetically cooled and 
trapped. Loss rates are obtained by exponential fitting to the maxima of the 
HD + resonances in the secular excitation spectrum (solid line in Fig. 12(b)). In 
this way, a 0.01 s _1 background loss rate of HD + ions from the trap is obtained 
when both the IR and UV lasers are turned off. This loss is due to chemical re- 
actions between HD + ions and background gases. The observed background loss 
rate is fitted well by a single exponential decay, which rules out strong nonlinear 
dependence of the Be + fluorescence during secular excitation on the number of 
HD + ions. 

The spectroscopic signal used to produce the spectra in Fig. 13 is the molec- 
ular ion dissociation probability, obtained as the relative change of the heights of 
the HD + secular resonances in the Be + fluorescence before and after the REMPD 
excitation (Fig. 12(d)). For each transition, the HD + dissociation probability was 
measured as a function of the frequency of the IR laser, in steps of 15 MHz. Each 
data point was obtained by averaging over several individual measurements of 
the HD + dissociation probability occurring over ~5 s. Each data point requires 
a new loading of HD + ions in the Be + crystal. For all measurements, comparable 
HD + ion numbers were used, as deduced from the size of the crystal core after 
loading. However, during each HD + loading cycle a small fraction of the Be + is 
lost from the trap, due to chemical reactions with neutral HD gas [63] . The same 
Be + ion crystal can be used for up to 40 HD+ loadings, sufficient for obtaining 
the spectra in Fig. 13. A typical spectrum is taken within 1-2 hours. 
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Fig. 13. Rovibrational transition spectra with partially resolved hyperfine splitting: 
(a) (v = 0,L = 2) -> («' = 4, L' = 1) at 1430 nm, (b) (v = 0, L = 2) -> («' = 4, 1/ = 3) 
at 1394 nm. The curves are fits to the data (•), where the theoretical stick spectra 
were broadened by «40 MHz. The theoretical spectrum exhibits a large number of 
very weak transitions, due to weak mixing of pure coupled angular momentum states. 
The ordinate values are the molecular ion dissociation probability for a 5 s irradiation 
of 0.65 W/cm 2 IR and 10 W/cm 2 UV light. The insets show typical error bars [69]. 

Detailed measurements for two transitions (v = 0, L = 2) — » (i/ = 4, L' = 1, 3) 
are shown in Figs. 13(a,b). Both spectra reveal a partly resolved hyperfine struc- 
ture, which can be compared with the prediction from an effective spin Hamil- 
tonian, written as H cS = bil p ■ S + ciI pz S z + b 2 I c i • S + c 2 IdzS z + 7S • J [72,73]. 
Here, I p , Id, and S denote the spin of the proton, deuteron, and electron, re- 
spectively; the subscript z indicates the projection on the internuclear axis. The 
hyperfine coefficients b\, 62, c\, C2, and 7 have been recently calculated to high 
accuracy [43], see Sec. 2.4. The hyperfine level energies and eigenfunctions are 
found by diagonalization of the matrix representation of -ff e ff in a suitable an- 
gular momentum coupling scheme. Terms arising from the nuclear spin-rotation 
and deuteron quadrupole interactions are neglected as they contribute -Cl MHz 
to the hyperfine level energies [43] . The results of the diagonalization were sub- 
sequently used to calculate line strengths (Eq. 11) of the individual hyperfine 
components within a given rovibrational transition, leading to " stick spectra" , 
as shown in Fig. 13. Inhomogeneous broadening of the spectra may be accounted 
for by convolving each line with a Gaussian lineshape of a given width. 

The broadened stick spectra are fitted to the experimental spectra using the 
lincwidth, the vertical scale and the frequency offset as fit parameters 
(Fig. 13). The frequency offset corresponds to the deperturbed ro- vibrational 
transition frequency, which is thus determined to within the accuracy of the 
wavelength calibration of the IR laser (40 MHz) and the fit uncertainty (3 MHz). 
The measured deperturbed ro- vibrational transition frequency is in good agree- 
ment with the ab initio results from [21], see Fig. 13. The partly resolved hyper- 
fine structure in the measured spectra agrees well with the theoretical results 
obtained from [72,43]. We find both theoretically and experimentally that the 
hyperfine structure for other transitions in the P and R branches is similar to 
that in Figs. 13(a,b). 
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We observe a typical line broadening of 40 MHz, which corresponds to fcs x 
(0.2 K) of energy in the axial motion. The kinetic energy in the secular motion 
(as inferred from molecular dynamics simulations) of the HD + ions can give 
rise to broadening of about 10 MHz only [61]. Saturation broadening also does 
not play a significant role, as confirmed by comparing spectra taken at different 
IR and UV intensities. Using the polarization-dependent 313 nm fluorescence 
of the Be + ions as a magnetic field probe, the magnetic field (which is along 
the direction of propagation of the 313 nm laser beam) has been adjusted and 
verified to be 50 mT and to vary by no more than 40 mT over the extent of 
the crystal, which implies Zeeman broadening <1 MHz. This leaves Doppler 
broadening due to micromotion as the most probable cause for the observed line 
broadening. This micromotion could arise from phase shifts in the rf potentials 
applied to the various electrodes, and from coupling between axial (IR laser 
beam direction) and radial ion motion. For our trap, in which the HD + ions are 
located at least 10 ^.m away from the trap axis, the (radial) micromotion energy 
exceeds fcs(0.5 K). 

The results described are of significance in several respects. They demon- 
strate, for the first time, the possibility of high-resolution spectroscopy of small, 
trapped molecular ion samples, sympathetically cooled well into the millikelvin 
range. We have achieved a spectral resolution 10 times higher than with any 
previous cold molecular ion method, and the same enhancement was obtained 
for the excitation rate. The observed population dynamics demonstrated the 
weakness of collisions. The methods used for trapping, cooling and detection are 
quite general, and are applicable to a host of other molecular ion species. This 
includes other ions of astrophysical and cosmological interest such as H^ and its 
isotopomers, which have been trapped in our setup [61,62]. Also, the spectral res- 
olution achieved here may be further improved: for instance, first-order Doppler 
broadening may be circumvented by use of a tightly confining trap which holds 
the ions in the Lamb-Dicke regime, or by two-photon spectroscopy. Furthermore, 
the presence of the atomic coolant ions offers an in situ tool to detect possible 
perturbing fields. 

5 Conclusion and Outlook 

In summary, the development of high-accuracy laser spectroscopy of trapped 
MHIs has made significant progress. On the theory side, the energies have been 
calculated with a relative accuracy of the order of 1 ppb. Detailed predictions of 
the line strengths of one- and two-photon transitions have been given, which are 
important guides for the experiments. Certain systematic shifts (dc and ac Stark 
shifts [15,22]) have also been calculated, but are not described here. On the ex- 
perimental side, several important techniques have been demonstrated: cooling 
of MHIs to tens of mK, vibrational-state selective photodissociation, one-photon 
vibrational spectroscopy with spectral resolution at the level of 2-10 -7 , rotational 
population measurement, in-situ ion detection, tunable, high-power, continuous- 
wave narrow-linewidth laser for two-photon spectroscopy. Based on the present 
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results, it is expected that the two-photon spectroscopy experiment will ulti- 
mately allow a spectral resolution at the level of 3 • 10~ 10 , while the one-photon 
1.4 /xm HD + spectroscopy in the current apparatus will be limited by Dopplcr 
broadening to several parts in 10 8 . One-photon spectroscopy of HD+ vibrational 
transitions having longer wavelength or the use of a trap with stronger con- 
finement should allow reaching the Lamb-Dicke regime, with a strong increase 
in spectral resolution. As described above, two-photon spectroscopy is another 
alternative. 

For both ion species, the investigation of systematic shifts will become an 
important task. It is expected that Zeeman shifts and Stark shifts can be reduced 
or measured to a level below one part in 10 10 in a cold ion ensemble. This 
should enable comparisons of experimental and theoretical transition frequencies 
at levels below 1 ppb, and, in the longer term, the development of a novel 
approach to the measurement of mass ratios of electron and hydrogen isotopes. 

In the future, it may become attractive to use the method of quantum-logic- 
enabled spectroscopy [74,75]. Some of the experimental limitations (broad state 
population distribution, need for destructive detection of molecular excitation, 
systematic effects) encountered with the approaches described here could be 
substantially alleviated. 
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